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Abstract 

We study the critical properties of the monopole-percolation transition in U(l) 
lattice gauge theory coupled to scalars at infinite (f3 = 0) gauge coupling. We 
find strong scaling corrections in the critical exponents that must be considered by 
means of an infinite-volume extrapolation. After the extrapolation, our results are 
as precise as the obtained for the four dimensional site-percolation and, contrary to 
previously stated, fully compatible with them. 



Keywords: Lattice. Monte Carlo. Percolation. Critical exponents. Phase transitions. 

Finite-size scaling. QED. 

PACS: 11.15.Ha 12.20.-m 14.80.Hv 



1 Introduction 

The richness of the phase diagram of models containing charged scalars and/or fermions 
and abelian gauge fields has produced the hope of finding a non-trivial critical point in 
four dimensions. For the compact formulation of lattice QED, the different phases of these 
models (Confining, Higgs, Coulomb), can be characterized in terms of their topological 
content and this relation can be numerically explored |2| (see Ref. for a detailed 
exposition). The link between topology and phase-diagram is strong to the point that a 
monopole-percolation second-order phase transition can be found beyond the end-point of 
the Higgs-Coulomb phase transition line (see Ref. [[| for a study of the phase-diagram). 
In Ref. ||, it has been conjectured that the monopole-percolation phase transition can 
produce a chiral phase transition that, being driven by the monopole percolation, would 
present the same critical exponents. 
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This conjecture has been put to test in Ref. ||, where the position of both the 
monopole-percolation and chiral critical lines have been located, by means of a Monte 
Carlo simulation in the quenched approximation. The chiral critical line is very close 
to the monopole-percolation one, but they can be clearly resolved in the limit of infi- 
nite (P = 0) gauge coupling. A clear-cut test of the scenario proposed in Ref. would 
be an accurate measure of both the chiral and monopole-percolation critical exponents. 
The critical exponents for monopole-percolation have been measured in Refs. || ||. The 
critical exponents displayed a very mild variation along the critical line, consistent with 
a single Universality Class, although significantly different from the site-percolation 0. 
However, in Refs. || [f| the scaling corrections are not considered in the analysis. 

In this letter, we report the results of a Monte Carlo calculation of compact scalar QED 
in the strong coupling limit for the gauge field (J3 = 0), where the model is integrable. In 
this way, we are able to directly generate independent configurations obtaining accurate 
measures in large lattices. 

In our study, we shall use a Finite-size Scaling (FSS) method based on the comparison 
of measures taken in two lattices, at the coupling value for which a renormalization group 
invariant (namely the correlation length in units of the lattice size) takes the same value 
in both lattices 0, ||, ||. In this way, by considering the scaling of other dimensionless 
quantities, we obtain direct information on corrections to scaling. In the present case 
the scaling corrections are notoriously difficult to deal with. The usual strategy of just 
considering the leading scaling corrections would only work in extremely large lattices, 
and we are compelled to use sub-leading corrections in the infinite volume extrapolation. 
After this extrapolation, we conclude that monopole-percolation belongs to the same 
Universality Class that site-percolation. 

2 The model 

The action for the (compact) U(l)-gauge model on the lattice, coupled to unit modulus 
scalars can be written as 

S = -(3J2 ReC/^(r)-«^Re$t( r )c/ M ( r )$( r + /A ) , (i) 

r,fi<u r,fi 

where r is the four-dimensional lattice site, /i and v run over the four spatial directions, \x 
is the vector joining neighbours along the \i direction, U^ u (r) is the elementary plaquette, 
Ufj,(r) the gauge variable, and $(r) the scalar field. The lattice volume is V = L A and 
periodic boundary conditions are imposed. Notice that we use the normalization of the 
parameter k as in Ref. |§] which is twice that used in Ref. ||. 

For (3 = and after a gauge transformation to eliminate the $ fields, the action 
simplifies to 

S = -K^ReC/» . (2) 

r,fj, 

The generation of independent Monte Carlo configurations for this action is straightfor- 
ward, as the link variables are dynamically independent. 
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To study the monopoles in the lattice, let us write U^(r) = e 19 ^ and define 

W) = W + 6 ^ r + M) - + ") - = WO + 2*N^(r) , (3) 

where 9^ is taken in the interval (— tt, ty] and is an integer. 
We obtain the monopole current in the dual lattice as @ 
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-e^ pa ^N pa (r + fi) , (4) 
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where A + is the forward difference operator in the lattice. Each component of the current 
m M is an integer which lives in the link of the dual lattice leaving r in the fi direction. 
Clusters are defined as sets of sites of the dual lattice connected through links with nonzero 
monopole current (occupied links). 

The observables that we measure for every gauge configuration are the link energy 
and magnetization-like quantities that can be expressed in terms of the cluster-size dis- 
tribution: 

E = ^Ref^(r), 

M 2 = Y^nl, (5) 

c 

M 4 = 3Mf - 2 n t , 

c 

M max = maxn c , 

c 

where n c is the number of sites of the c— nth cluster. E is used to compute k derivatives 



and to extrapolate the measures taken at k to neighbouring coupling values fL0 |. The 
definition of M 2 (M 4 ) can be understood by putting in the occupied sites Ising spins 
at zero temperature (spins on the same cluster have the same sign), taking the second 
(fourth) power of the magnetization, and averaging over the signs of the clusters 

The same construction allows for a sensible measure of the correlation length in a 
finite lattice. We first measure the Fourier transform of the clusters 



n c (k) = J2 eik ' r ( 6 ) 

r£c 

at minimal momentum, from which we obtain 

Ekwi 2 ) > ( ? ) 




and then use the following definition [[□]] 



(M 2 )/F-l\ 1/2 



4sin 2 (7r/L) 
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We have used two definitions of the unconnected susceptibility 



Xi = (M 2 )/V, 

X2 = {MlJ/V. W 

Notice that the monopole density M\/V is not critical at the transition. So, we can define 
also the susceptibilities dividing M 2 or M max by Mi, as done in Refs. || |(|. This should 
only modify the corrections to scaling (in general, we find that they increase slightly). 

The previous numerical studies have been mainly based on measures of connected 
susceptibilities, specifically 

Ti = ((M 2 -Ml ax )/M 1 ) 

% = (((Mnax/Mx) 2 ) — (M mSLX /Mi) 2 )V . [W) 

These definitions make sense at both sides of the transition and present a peak near it. 
However we shall see that they present strong corrections to the scaling, becoming less 
appropriate than the unconnected ones for a FSS study. 

It is also very useful to measure quantities that keep bounded at the critical point, 
but whose k derivatives diverge. Some examples are the correlation length in units of the 
lattice size and the Binder parameters: 

1 ( (M4} 



2 V~ (M2) 2 
M 

im 

(M 2 ) 

\ max / 



z 

b W fin 



3 The numerical method 

We have simulated in symmetric lattices of linear sizes L = 6, 8, 12, 16, 24, 32, and 48. We 
have generated 10 6 independent configurations for each lattice. The statistical analysis 
have been done with 1000 bins of data, for an accurate error determination. To speed 
up the computations, we have used the U(l) subgroup ^5535. We have checked that the 
U(l)-discretization effects are negligible comparing the results with those using Z255. 

In order to measure critical exponents, we use a FSS method. Specifically, we use the 
quotients method used in Refs. 0, f| ||. Given an observable O that diverges as t~ x °, 
t being the reduced temperature (k — k c )/k c , the FSS ansatz predicts that for a finite 
lattice of size L, in the critical region 

(0(L,t)) = L*°'"[F (Z(L,t)/L)+0(L-»)] , (12) 

where Fo is a smooth scaling function and u is the universal exponent associated to the 
leading corrections to scaling. 
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To eliminate the unknown function Fo, one can compute the quotient Qo of the mean 
value of the observable in two different lattices, at the coupling value where the correlation 
lengths in units of the lattice size is the same (there is a crossing): 

«°i*-= fra=^" +0 < £ ->- (13) 

lattice sizes being sL and L respectively. The 0(L~ W ) terms include all powers of as 
well as the corresponding series produced by sub-leading irrelevant operators || [12|, [13|] . 
Other corrections are generated by the non singular part of the free energy (analytical 
corrections). For most observables the analytical corrections are 0(L~ 7//l/ ). Let us remark 
that Qo and are statistically correlated, which allows for an important (statistical) 
error reduction in Eq. ([13]). In fact, a modification of Eq. (|13|) can be used to study loga- 
rithmic corrections to mean-field behaviour with rather high accuracy ]14j . Let us finally 
remark that a similar two-lattices matching method has also being extremely successful 



in lattice QCD studies [15 



It will be also useful to recall the shift of the finite-lattice critical point from the true 



critical point 16 



ros S(LjS)(x i_^ L -a,-l/, ) (M) 

where only the leading scaling corrections have been kept. With a fixed value of s, 
t cross (L, s) oc J_,- W - 1 I V ^ to be compared with the shift of the peaks of the connected suscep- 
tibility, that goes as 

t peak K L -l/u Thege peakg 

can be measured in a single lattice, but 
we loose a factor of . That is why the quotient method suffers from smaller scaling 
corrections. Eq. ([14]) also applies to the crossing of the Binder cumulants Bi. 



4 Results 

In table [l] we show our results for the critical exponents 77 and v as obtained from Eq. fllTf). 
The used operators have been xi ( x xi = 1 jv = 1 — r\) and <9£/cta (xq k ^ = 1 + 1/u). In 



v 77 Z/s-p ?7s-p 

6 0.6801(34) -0.2336(7) — — 

8 0.6829(36) -0.1713(6) 0.689(3) -0.0687(7) 

12 0.680(5) -0.1256(6) 0.687(3) -0.0775(7) 

16 0.681(6) -0.1095(7) 0.688(4) -0.0825(6) 

24 0.689(6) -0.0986(6) 0.691(5) -0.0868(8) 



Table 1: Critical exponents obtained from Eq . (|T3|) with lattice pairs (L,2L). We have 
used as operators <9£/<9k (xq k ^ = l + l/u), and xi ( x xi = ll v = 2 — 77). The last two 
rows display the corresponding results for site-percolation (obtained with exactly the same 
operators) reported in |7]]. 
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all cases, we consider only the ratio s = 2. These exponents can be directly compared 
with the results for the four dimensional site-percolation 0. The trend for exponent 
v in both cases is rather similar: the scaling corrections are significantly smaller than 
the statistical errors, so that the results seems stable with growing L. To this accuracy, 
both measures are compatible. On the other hand, for the anomalous dimension, rj, 
both systems present significant scaling corrections, although stronger in the monopole- 
percolation case. Therefore, the scaling corrections must be dealt with before comparison 
for 77 can be attempted. 

In previous studies @, H, we have used Eq. ( |14|) to obtain an estimate of uj that 
allows to perform an infinite-volume extrapolation. However, in this case the higher- 
order scaling corrections are so large that they need to be considered. This can be seen 
in figure |2], where we show the finite-lattice critical point as obtained by the crossing of 
£/L, and of Bi. In the £/L case, the behaviour is not even monotonous with increasing 
lattice size. Therefore, in this problem we need to go beyond Eq. (14"). 



Thus, we have considered the three Binder cumulants B { whose quotients should be- 
have as 

Qb^s = 1 + AL-" + B t L- 2 " + CiL'"' + D.L-^ + . . . , (15) 

where uj' stands for a sub- leading irrelevant exponent. 

Our numerical results for these quotients are plotted in figure [L| Although the be- 
haviour with growing L is not monotonous, we have first considered the parametriza- 
tion 1 + AiL~^ , which should be adequate for large L. The quality of the fit is rather 
poor unless discarding all but the two largest pairs: uj = 0.81(2) (x 2 /dof = 159/8); 
to = 0.99(5) (x 2 /dof = 19.9/5); uj = 1.03(8) (x 2 /dof = 5.2/2), for L > 8, L > 12 and 
L > 16 respectively ("dof" is the number of degrees of freedom in the fit). We see that 
higher-order corrections need to be included. As uj seems slightly larger than one and, 
from table [I] it is clear that j/u will be close to two, L~ 2u and L~~'/ u will be of the 
same order. Therefore, a quadratic fit in L~ u should be a good parametrization provided 
that there are not sub- leading irrelevant operators in the intermediate range. In fact, the 
quadratic fit shown in figure p] (discarding the smaller lattice, L = 6) is quite reasonable: 

uj = 1.23(8), x 2 /dof = 5.4/5 . (16) 

Notice that the data used in the fits are strongly correlated and the consideration of the 
full covariance matrix is mandatory. 

A cross-check of the determination Eq. (|16D can be done fitting to the functional form 
1 + AiL-" + CiL^I\ fixing 7/1/ = 2.09 (this is the value found in Ref. [§, for site- 
percolation and it is also quite close to the values encountered in table [I]). This fit yields 
a compatible value with a rather increased error: uj = 1.36(14). 

Therefore, our estimate seems consistent although one may prefer to double the error 
in Eq. (Jl6|) , as an estimation of systematic errors, to be in the safe side. 

A different cross-check is the infinite- volume extrapolation for k c , using Eq. (Tl) and 



the crossing points for £/L, Bi,B 2 and B 3 . In figure we plot the cumulant crossing 
points as a function of L~^~ x l v . As the behaviour is not monotonous we have tried a 
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Figure 1: Quotients of cumulants B\, Bi and B3 as a function of . We use for u the 
value obtained in the fit, thus the curvature is related with higher-order corrections to 
scaling. 



fit to k c + AiL-"- 1 /" + BiL- 2w - l / u . From table [I] and Eq. (|IJ) we find that v can be 
determined much more accurately than lu, so its value can be safely fixed to v = 0.69 in 
the fit. Fixing also u to (|16|) we obtain an acceptable fit only if L > 12 (x 2 /dof = 1.1/3). 
We obtain 

k c = 2.698736(34)(11) . (17) 

Through out the paper the second error will denote the error induced by the uncertainty 
in u. Therefore, if one chooses to double the error in u, this second error needs to be 
doubled too. As the value of k c will be by far our most precise result, it is important to 
check that the discretization effect of using Z G55S5 instead of U(l) is negligible. In order to 
do so, we have repeated the measure of the crossing point for the pair (12, 24) using Z255 
with the same statistics, obtaining a compatible value within errors (one per million). 

The exponent r\ can be obtained from the susceptibilities \\ an d X2 (which scale as 
jjilv = L 2-^ and the ma g ne ti za tion (M max ) (which scales as L~ p/u+4 = L 3 "^ 2 ). In 
figure H we show the three rj determinations. It is clear from the plot that we cannot 
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Figure 2: Crossing points of the cumulants, as functions of L w x l v (here we have set 
u = 1.23, v = 0.69). 



keep just the leading scaling corrections. In this figure we show a joint fit of the three 
quantities, quadratic in L~^, with L > 8. The extrapolated value is 



n 



-0.0876(22)(6), x 2 /dof = 1.2/5 . 



This might be compared with the result obtained by assuming that the scaling corrections 
are A { L^ + QL^ F (fixing again 7/1/ = 2.09): 



n 



-0.0902(25)(1), x 2 /dof = 2.4/5 . 



(19) 



We see that by following our recipe of doubling the uj induced error, both extrapolations 
are not covered by the second error bar. We thus take the difference to estimate the 
systematic error involved in the infinite-volume extrapolation generated by higher-order 
terms. The final value that includes both kind of errors is 



n 



-0.089(4) . 



(20) 



For computing the v exponent we measure the quotient corresponding to the operator 
d^/dn, which scales as L 1+1 ^ u . From table [l] we see that the infinite- volume extrapolation 
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Figure 3: The exponent rj as a function of for several magnetization or susceptibility 
operators. The dashed lines are quadratic fits in the variable L~ w (u = 1.23). Notice 
that the connected susceptibilities present stronger scaling corrections. We also plot the 
results of the site-percolation taken from Ref. 0. 

is not so crucial in this case. However, we cannot just average the different determinations 
(which is basically what a log- log fit does), as the statistical error in the mean can decrease 
enough to uncover the scaling corrections. In order to obtain a safe error estimate, we 
perform a fit linear in . The difference between the extrapolation with L > 6 and with 
L > 8 is ten times smaller than the error. Taking the error from the latter we obtain: 

v = 0.685(6), xVdof = 1.0/3 . (21) 

In this case, the cu-induced error is twenty times smaller than the statistical error. 

We summarize the infinite- volume extrapolation of critical exponents and of the critical 
coupling in table |[ 
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5 Conclusions 



We have studied with Finite-size Scaling techniques the monopole percolation transition 
of compact QED coupled to scalars in the strong-coupling limit ((3 = 0). We have shown 
that state of the art techniques for measuring critical exponents in spin models can be 
successfully applied to this problem. The approach relies in comparison of measures taken 
in two different lattices when a matching condition is fulfilled (see Eq. (|13|)). The efficiency 
of the method is greatly enhanced by the availability of a re- weighting method ||10|| , and 
of an easily measurable Renormalization Group invariant (the correlation length in units 
of the lattice size [PH ). 

With the achieved accuracy in individual measures, the scaling corrections for the 
reachable lattice sizes are big enough to require the consideration of sub-leading scaling- 
corrections in the infinite-volume extrapolation. After this extrapolation, it is found 
(within errors) that at (3 = monopole-percolation and site-percolation belong to the 
same Universality Class. The discrepancy with previous calculations is explained by the 
presence of strong corrections to scaling. 

The present study can be extended to the monopole-percolation critical line at non- 
zero (3 coupling [[| although the number of independent configurations that could be 
generated would be quite smaller. Another interesting matter is the comparison with 
the chiral critical behaviour. Our method would be useful in this respect only if it is 
found an analogous of the correlation length in a finite lattice, that made sense at (chiral) 
criticality. 
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Source 




V 


V 


UJ 


This work 


2.69874(6) 


0.685(6) 


-0.089(4) 


1.23(16) 


Franzki, Kogut, Lombardo fU 


2.6938(8) 


0.61(4) 


-0.28(2) 




Site-percolation R] 




0.689(10) 


-0.094(3) 


1.13(10) 



Table 2: Summary of the values obtained for the critical coupling including statistical 
and our estimation of systematic errors. In the second row we show the results of Ref. || 
for the same system. Finally (third row) we recall the results of Ref. |7j for the site- 
percolation. 
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